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Abstract 



In this paper, we introduce a method of imposing asymmetric conditions on the velocity 
vector with respect to independent variables and a method of moving frame for solving the three 
dimensional Navier-Stokes equations. Seven families of non-steady rotating asymmetric solutions 
with various parameters are obtained. In particular, one family of solutions blow up at any point 
on a moving plane with a line deleted, which may be used to study turbulence. Using Fourier 
expansion and two families of our solutions, one can obtain discontinuous solutions that may be 
useful in study of shock waves. Another family of solutions are partially cylindrical invariant, 
contain two parameter functions of t and structurally depend on two arbitrary polynomials, 
which may be used to describe incompressible fluid in a nozzle. Most of our solutions are 
globally analytic with respect to spacial variables. 

1 Introduction 

The most fundamental differential equations in the motion of incompressible viscous fluid 
are Navier-Stokes equations: 

Ut + UU^ + VUy + WU^ + = Viu^x + Uyy + ^^22), (1-1) 

Vt + UV^ + VVy + WV^ + -^Py = P{V^^ + Vyy + f^^), (l.2) 

wt + uw^ + vwy + ww^ + -^p^ = iy{w^^ + Wyy + w^^), (1.3) 
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Ux + Vy + Wz^O, (1.4) 

where {u, v, w) stands for the velocity vector of the fluid, p stands for the pressure of the 
fluid, p is the density constant and u is the coefficient constant of the kinematic viscosity. 

The Lie point symmetries of the two-dimensional special case of the above equations 
{uz — Vz — w — 0) were obtained by Pukhnachev [PI] and Buchnev [Ba]. Moreover, 
certain group- invariant solutions were found in the works of Pukhnachev [Pvl], Kochin- 
Kibel'-Roze [KKR] and Bytev [Bvl], [Bv2]. Puthermore, Gryn [G] obtained certain exact 
solution describing flows between porous walls in the presence of injection and suction at 
identical rates, and Polyanin [Pa] used the method of generalized separation of variables 
to find certain exact solutions. 

Assuming nullity of certain components of the tensor of momentum flow density. Lan- 
dau [LI] found a exact solution of Navier-Stokes equations (1.1)-(1.4), which describes 
axially symmetrical jet discharging from a thin pipe into unbounded space. The Lie point 
symmetries of the above three-dimensional equations were obtained by Buchnev [Ba] and 
Pukhnachev [Pv2] . Moreover, Kapitanskii [K] found certain cyhndrical invariant solutions 
of the equations and Yakimov [Y] obtained exact solutions with a singularity of the type 
of a vortex filament situated on a half hue. Shen [SI, S2] rewrote Navier-Stokes equations 
in terms of complex variables and found certain exact solutions. Brutyan and Karapivskii 
[BK] got exact solutions describing the evolution of a vortex structure in a generahzed 
shear flow. Furthermore, Leipnik [Lr] obtained exact solutions by recursive series of dif- 
fusive quotients, and Vyskrebtsov [V] studied self-similar solutions for an axisymmetric 
flow of a viscous incompressible flow. 

From algebraic point of view, it seems to us that there are not enough exact solu- 
tions that fully reflect the fundamental natures of Navier-Stokes equations. In this paper, 
we introduce a method of imposing asymmetric conditions on the velocity vector with 
respect to independent variables and obtain two families of non-steady asymmetric so- 
lutions with rotation. One of the families contains two arbitrary parameter functions of 
t and an arbitrary number of parameter constants. Using Fourier expansion and this 
family of solutions, one can obtain discontinuous sohitions that may be useful in study of 
shock waves. Another family of solutions are partially cylindrical invariant, contain two 
parameter functions of t and structurally depend on two arbitrary polynomials, which 
may be used to describe incompressible fluid in a nozzle. In order to better reflect the 
rotating nature of flow, we also give a method of moving frame and flnd flve families of 
non-steady rotating solutions with various parameters. In particular, one family of solu- 
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tions blow up at any point on a moving plane with a line deleted, which may be used to 
study turbulence. Another family can be used to obtain discontinuous rotating solutions. 
Most of our solutions are globally analytic with respect to spacial variables. Below we 
give a more detailed introduction. 

The equations (1.1)-(1.4) are invariant under orthogonal transformations {Ta \ A G 
0{n,R)} with 







Moreover, they arc invariant under the time translation t 
scahng with 7^ 6 e M: 



A, Ta{p)=p. (1.5) 
> t + a with a e M and the 



Tb{u) = b~\{bH, bx, by, bz), Tb{v) = b~^v{bH, bx, by, bz), 
Tb{w) — b~^w{bH, bx, by, bz), Th{p) = b~^p{b^t, bx, by, bz). 



(1.6) 
(1.7) 



The most interesting symmetries of Navier-Stokes equations are the following time-dependent 
translations: 

Tia{u) = u{t,x + a,y,z) - a', Tia{v) = v{t,x + a,y, z), (1.8) 

Tia(w) — w(t,x + a,y, z), Tia(p) — p(t,x + a, z) + pa''x (1-9) 
and its permutations on {u,x), {v,y), {w,z), and 

T2a{v)^v, T2a{w)^w, T^^ip) ^ p + a, (1.10) 

where a is an arbitrary function of t. The above transformations transform solutions 
of Navier-Stokes equations into their solutions. Our goal in this paper is to find exact 
solutions of Navier-Stokes equations modulo the above symmetries. In other words, the 
above symmetries will be used to simplify our ansatzes for exact solutions and related 
arguments. 

For convenience, we always assume that all the involved partial derivatives of related 
functions always exist and we can change orders of taking partial derivatives. In fluid 
dynamics, rotation-free solutions of Navier-Stokes equations, namely. 



Vx = 0, 



Wy = 0, 



Uz = 0, 



(1.11) 



are not so interesting. Prom pure mathematical point of view, a rotation-free solution is 
equivalent to a time-dependent harmonic function f{t,x,y,x) (i.e., f^x + fyy + fzz = 0), 
where 

u-^ fx, v = fy, w = fz. (1.12) 



Practically, steady solutions (or time-independent) are not very important. In general, it 
is difficult to find exact non-steady rotating solutions for Navier-Stokes equations (1.1)- 
(1.4) due to their nonlinearity. 

Using certain finite-dimensional stable range of the nonlinear term, we found in [XI] a 
family of exact solutions with seven parameter functions for the equation of nonstationary 
transonic gas flows found by Lin, Reisner and Tsien [LRT] , which blow up on a moving line. 
These solutions may reflect partial phenomena of gust. In [X2], we use various ansatzes 
with undermined functions and the technique of moving frame to find basic solutions 
modulo the Lie point symmetries with parameter functions for the classical non-steady 
boundary layer problems. These two works motivated us to solve Navier-Stokes equations 
by algebraic methods. 

Our first idea is to impose suitable asymmetric conditions on the velocity vector with 
respect to independent variables. For instance, assuming 

« = 7(t)x + #(t,x2 + ?/2), v^-i{t)y-x(t>{t,x''+y''), w ^ il;{t,x'' + y^) -2j{t)z, (1.13) 

we obtain the following solution of Navier-Stokes equations (see Theorem 2.4): 

a' Py ^ {adtYi^) f + y^ 



u = — X 



2a x'^ -\- y"^ ^ -\- l)\a \ Aua 



V-^c^^tM:^ (114) 



1=0 

, _ _ , V ^"^^^^(^^ (^^] ^ (1 15) 

2a^ x^ + y^ ^il{i + l)\a \ 4ua J ' ^ ' 

w^ay — — — z, (1-16) 

p((a'f - 2aa")(x^ + y^) p(aa" - 2(a'f)z^ ^, y 
p = ^-^^ ^ + ^1^ :r—^ h/?'arctan^ 

where a, (5 are any functions in t and $5, ip are arbitrary polynomials in t. The above 
solution can be used to describe incompressible fiuid in a nozzle. The polynomials and 
</? can be replaced by the other functions as long as the related power series converge. 

As we emphasized earlier, people are interested in solutions that are not rotation free. 
To better capture the rotating nature of fiuid, we introduce the following moving frames: 

A' = xcosQ! -|- (y cos/3 -|- 2;sin/3) sina, 3^ = — xsin a -|- (ycos/3 -|- zsin/?) cosa, (1-18) 

Z = —ysmP + z cos j3, U = ucosa+ {y cos P + wsmj3)saia, (1.19) 



V — —usina + {vcos (3 + wsin P) cosa, W — —vsinP + wcosl3, (1-20) 

where a and (3 are functions in t. Here we exclude the translation components because 
Navier-Stokes equations are invariant under the transformations of the type Tia in (1.8) 
and (1.9), and we want to consider sohitions modulo these transformations. With respect 
to the above rotating frames, Navier-Stokes equations change to more complicated system 
of partial differential equations. Imposing asymmetric conditions on the moving frames, 
we find another five families of non-steady rotating solutions with various parameters. For 
instance, we have the following solution of Navier-Stokes equations (see Theorem 3.3): 

u = + 6uyX-'^^ sin a - A" cos a) — a'{X sina + y cosa), (1-21) 

V — — (— — \- 6i'yX~'^] {X sina + y cosa) cos P + a'{X cosa — ysina) cos P 
\2a' J 

-P'ZcosP+ (^l3'Xsma + I3'y cos a- sin P, (1.22) 



w 



— ( ^ — h GuyX ^ ) (Xsina + y cos a) sin P + a'(X cos a — 3^ sin a) sin P 
\2a' ) 

-P'ZsmP+ (^Z-P'X sin a — P'y cos cos P, (1-23) 

^ i2a'a'" + 4{a')^ - Sja")^) {X^ + y^) ?,{p'f{X'' sin^ a + y cos^ a) 

P - %{a'Y 2 

+12iy(a'yX-^ - uX-'^) + {p" - 0^)Z(X sin a + :V cos a) 

.^/N2..., (a'a'" + (a'y(p'y -2(a"f)Z\ , 

-ZiP'fXy sin a cos a + ^ ^ — '-^ — }. 1.24) 

2\a!Y 

The above solution blows up at any point on the following rotating plane with a line 
deleted: 

{(a;, z) e I a; cos a + ?/ sin a cos P + z sin a sinP — 0, 

—xsina + ycosa cos/3 -|- 2; cos a sin/3 7^ 0}. (1-25) 

This type of solutions may be applied in studying turbulence. Since all of our solutions 
in this paper only involve elementary functions and integrations, they may be applied to 

engineering problems with the help of computer, although they appear sophisticated in 
format. They can also be used to solve certain initial value problems for Navier-Stokes 
equations because they contain parameter functions. 
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As we all know that in general, it is impossible to find all the solutions of nonlin- 
ear partial differential equations both analytically and algebraically. In our arguments 
throughout this paper, we always search for reasonable sufficient conditions of obtaining 
exact solutions. For instance, we treat nonzero functions like "nonzero constants" for this 
purpose because our approaches in this paper are completely algebraic. Of course, one 
can use our methods in this paper to get more solutions, in particular, by considering the 
support and discontinuity of the related functions. We want to remind the reader that 
we always put arguments (proofs) before our theorems (conclusions) due to our purpose 
of finding exact solutions. 

The paper is organized as follows. Section 2 is devoted to our asymmetric approaches. 
We present the general settings for the moving-frame approach in Section 3 and find two 
families of exact solutions. In Section 4, we use the moving frames and certain ansatzes 
involving given irrational functions to find another three families of exact solutions. 

2 Asymmetric Approaches 

In this section, we will solve incompressible Navier-Stokes equations (1.1)- (1.4) by impos- 
ing asymmetric assumptions on u, v and w. 
For convenience of computation, we denote 

^l^Ut + UUx + VUy -\- WUz - u{Uxx + Uyy + u ^ z) , (2.1) 
^2^Vt + UV^ + VVy + WVz - y{Vxx + Vyy Vzz) , (2.2) 

^3 = wt + uw^ + vwy + wwz - y{wxx + Wyy + Wzz)- (2.3) 
Then Navier-Stokes equations become 

*i + -Px = 0, $2H--Py = 0, $3 + lp, = (2.4) 
P P P 

and Ux-\-Vy-\-Wz — 0. Our strategy is first to solve the following compatibility conditions: 

a,($i) = a,($2), 5,($i) = a,($3), 5,(^2) = ^,($3) (2.5) 

and then find p via (2.4). 

Let us first look for simplest non-steady solutions of Navier-Stokes equations (indeed, 
the corresponding Euler equations) that are not rotation free. This will help the reader 
to better understand our later approaches. Assume 

u = 'jix - aiy - a2Z, v = aix + 72^/ - a^z, w = + a^y + 73Z, (2.6) 
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where cuj and 7^ are functions in t such that 71 + 72 + 73 = 0. Then 

*i = (Ti + 7i - - OiD^ - (ct'i - «i73 + Oi2a3)y + {aias - ol^ + 0:272)2;, (2.7) 

^2 = (ai - q:i73 - Q;2a3)a; + (72 + 72 - «i " "3)?/ " (^Ja + «ia2 - ^371)^, (2-8) 

^3 = (a2 + OL\OL-A - OLll2)X + (Og - Q;iQ;2 - Q!37l)Z/ + (73 + 73 " "2 " "s)^" (2-9) 

Furthermore, 



5,(*i) = d,{^2) ^ 73 = — , (2.10) 



a' 



a,($i) = a,($3)^72 = — , (2.11) 

q;2 



a' 



a,($2) = a,($3)^7i = — • (2.12) 

Ct3 



Note 



«'i CKo a'q 



71 + 72 + 73 = \ \ = 0~ a.\a.2(y.'i = c (2.13) 

(X\ OL2 OL-i 



for some real constant. Moreover, 

$1 = (oig'ttg ^ — al — a\)x — a20ii,y + a-ia^z, (2.14) 

$2 = — Q;2Q;3a: + {0.20-2^ — a\ — a'^y — aia2Z, (2.15) 

$3 = a^a^x — aia2y + {ot^al^ — a\ — a^)z. (2.16) 

By (2.4), 

P\/ 2,2 /' -i\ 2 I / 2 I 2 /' -i\ 2 I / 2 I 2 /' -i\ 21 

P — -[(Q;^ + — dig )X + (CK]^ + — CKg )y + (CKg + CVg — CK]^ q;^ J 

+p(q;2Q;3X?/ — aiazxz + Q;iQ;2y^) (2.17) 

modufo the transformation in (1.10). 

Proposition 2.1. Let ai, a;2 o'^t^ ctz be functions in t such that ai«2<^3 = c for 
some real constant c. Then we have the following solution of Navier-Stokes equations 
(1.1)-(14): 

a^' a2 Oil 

u = — x — aiy — a2Z, v = aix H y — a^z, w = a2X + a^y H z (2.18) 

0-2 Oil 

and p is given in (2.17). 
Next we assume 

^^ = -^2/, w^i:{t,z), (2.19) 
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where /3 is a function in t, ijj is a, function of t, z and v is so written just for computational 
convenience by our earlier experience in [XI, X2]. According to (1.4), 

u = fit,y,z) + -ip,^ X 
for some function f of t,y,z. Then 



+ 



W~ ) — 2^0^ 'I'^t-i^Az + i^i^z 



$9 



Thus (2.5) is equivalent to the following equations: 



r 



/ (3" \ (3" 
ft + f i-^pi ~^^] ~ Q^^y-fy + ~ ^^fyy + 



0, 



r 



(3" 



= 



with T = dy, dz- 

Given 6, c e M and a function 7 of we set 



^0 = he^^+^i - ce"^^-'^^, ii = b sin( ^/^z - v^) , 
Then we have the following solution of (2.24): 
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with r = 0, 1. Write 



ii'^W 2y 
/ 



The equation (2.23) is implied by the following equation: 



To solve the above equation, we assume 



(2.20) 



(2.21) 
(2.22) 

(2.23) 
(2.24) 



(2.25) 
(2.26) 

(2.27) 
(2.28) 

(2.29) 



f^h{t,y)+g{t,y)Cr 



(2.30) 
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where h and g are functions in t, y. Then (2.29) is imphed by the following two equations: 

(2.31) 



h 

9t - ^y9y - ^9yy " ^ = 0, 



(2.32) 



We will solve the above system of partial differential equations according to the following 
two cases: 

Case 1. r = 0, 6c = or r = 1, 6 = 



In this case, we have the following solutions: 



i=i J PI 

n 

h^J2 c/2,.e'^(«t-''i,J/^+«2-^^ sin(62,.(2z.a2,,/? + ^y) + C2,s), 



(2.33) 

(2.34) 



s=l 



(2.35) 



where ai^j, Ci^j, di^j and 02,5, 62,8) C2,s, c?2,s are real constants. 
Case 2. r — 0, be ^ or r — 1, 6 7^ if r = 1. 
For convenience, we denote 

a — y^4:5o,rbc + 5\^r^. 
Then we have the following solution of the system (2.31) and (2.32): 

g = cosh (^a J d,,ie^^<^-'h^^+'^^-^y sin(6i,,(2z.ai,/5 + ^'y) + Ci,,) 

sinh \ a f '^^ ) ^ 

+ ^ ^ Y rf2,.e^^«t-^t)/^+«^.^V3^^ sin(62,.(2i/a2,./3 + ^/^y) + C2,.), (2.36) 



s=l 



h^a sinh (^a J c/i,ie^(«M-''M)^+"^'^^^ sin(6i,,(2i/ai,,/? + ^/p^y) + Ci,,) 

+ cosh (^J |jd2,.e^^4^-^t)/^+«^.^VFy3in(6,^^(2^a2,./?+ V^?/) + C2,.). (2.37) 



In any case, 

_ . 2P'P"'-{P"f 2ii"-{i'f V46c+^i,,6^ 



4(/3')^ 



4(7')2 2/?Y (/?'7')^ 



(2.38) 



Thus (2.4), (2.38), (3.39) and the first equation in (2.22) give 



{2{3'I3"' -Z{(3"Y)y^ I 2y/^-3(yO' , (-1)^(2/3^ + 5/?Y)C. - C' 
^ ^ 8(/3')2 8(y)2 ^ 2(/3')^(7')' 

px^ / {I3"Y - 2(3' P'" , (V)^ - 2y7'" (3"i' 5o,Abc + 5,,rb' 



(2.39) 



+ 



4(y)^ 



2/3Y (/3'y)2 



(2.40) 



modulo the transformation in (1.10). 

By (2.20), (2.27), (2.28), (2.33), (2.34), (2.36) and (2.37), we obtain: 

Theorem 2.2. Let (3, 7 he functions in t, and let b, c, Oi^j, bi^i, Ci^j, di^i and 02,5, 62, si C2,s, 
(^2,5 be real constants. Define CLnd C,r in (2.25) and (2.26). For r = 0, 1, we have the 
following solution of Navier-Stokes equations (1.1)-(1.4): 



w 



(^'VW 2y 



(2.41) 



p is given (2.^.0), and 

c 



u 



'= d,,e^^<^-'h^P+'^^'^^y sm{b,,{2va,,(3 + ^'y) + ci,) 
X sin(62,,(2z/a2,,/5 + + C2,.) (2.42) 



+ 1,2/3' + 2y /3Yr h 



if r — 0, be — 0, and 

^)' + ^(^'"*(/^) """""" 

m 



^ ^2/3' 2y /3'7' 



i=l 



^(-(/^)4-(/^))|:*.«-'- 



2 -6i „)/3+a2,.V3"!/ 



X sin (62,3 (2i^a2,s/9 + V^?/) + C2,s) 



if r = 0, be ^ or r = 1, b ^ 0, where a is defined in (2.35). 



(2.43) 



Remark 2.3. We can use Fourier expansion to solve the system (2.31) and (2.32) 
for g{P, y/P'y) and /i(/3, y/P'y) with given g{0, \f^y) and /i(0, \f^y). In this way, we can 
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obtain discontinuous solutions of Navier-Stokes equations (1.1)- (1.4), which may be useful 
in studying shock waves. 



Set 

w^x^ + (2.44) 

Consider 

cx! q! q! 

—x + y(j){t,w), V ^ —y - x(j){t,zu), w^'4'{t,zu) z, (2.45) 

where a is a function in t and 0, are functions in t^w. Then (2.1)-(2.3) give 

2q;q;'' — (a')^ , "V/ ,v . , ,v ,^ 

' = 4q,2 ^ + + a ^"^^^"^ ~ ~ ^yyiwcf))^^, (2.46) 

$2 = , o y - a;0t (tJ70)^ - + 4a;i/(ro0X^^, 2.47 

2(0;')^ — aa'' , a' , a' , , , , , , 

$3 = p ^ + -0* + — ^^7^^^ - 4i/(V'^j7 + w^^^). (2.48) 

Note that dy{^\) = dx{^2) becomes 

a' 

{w(p)^t H {.{^(P)u7 + '^{'^(p)mm) - 4z/((q70)^j,^ + w{w4>)^^^) = 0. (2.49) 

a 

Set 

= Q;(ca(/))^. (2.50) 

Then (2.49) becomes 

a' - 

0t H txJ0^jj - 4i/(0^ + w4>^^) = 0. (2-51) 

a 

So we have the solution 

^ («a,)-(S) / y 

j=0 ^ ^ 

for a polynomial in t. By (2.50)-(2.52), we have 

^ = +Z.^!(-n)!^(4^j (2-53) 

1=0 ^ ' 



for a function (3 of t. 
Note 

oo 

0, = ^'w-' + 5^ 



1=0 



i!(i + l)!a2 (z!)2!a2 j ' (^.54) 

a' , a' 7 a' ^ (Q;9t)*($5) f w y 

= = E ^7^^ (i;;;;) ■ (2-55) 

i=0 ^ 
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Thus 



Therefore, 



=1 ^ ' i=o ^ ■' 

bt + -{w(t>)^ - Aiy{wct>)^^ = P'zu-\ (2.57) 



q; 



and 



2aa"-(a'y I3'y ,^ , , 

On the other hand, Equations dz{Ri) = dx{R3) and dz^R^) = dy{R3) are imphed by 
the following differential equation: 

a' a' 

i/^t H ruip-^ - 4zy(V'ro + Tuip^-^) = (2.60) 

a a 

(cf. (2.48)). Similarly, we have the solution: 

s=0 ^ 

where (p is a polynomial in t. With this ip, 

$3 = z. 2.62 

By (2.4), (2.44), (2.45), (2.53), (2.58), (2.59), (2.61) and (2.62), we obtain: 

Theorem 2.4. Let a, {3 he any functions in t and let </? he polynomials in t. We 
have the following solution of Navier- Stokes equations (1.1)-(1.4): 

« = — x+ , , +y> \, —] , 2.63 

2a x^ + y^ '^«K^ + l)!aV 4z/a / ' ^ ^ 

Px ^{am^(x^+jl_\ 
2J x^ + y^ ^^z!(z + l)!a V 4z^« 



piia'f -2aa"){x^ + y^) piaa" - 2ia'f )z^ ^, y 

p = ^-^^ ^ + ^ y ' ' h 13' arctan ^ 

80;"' x 



+2"^" i\ii + iy.s\{s + mar j ■ ^2.66) 

Remark 2.5. The above solution can be used to describe incompressible fluid in a 
nozzle. The polynomials $5 and f can be replaced by the other functions as long as the 
related power series converge. 



12 



3 Moving-Frame Approach I 



In this section, we will present the general settings for the moving-frame approach and 
find two families of exact solutions. 

Let a, 13 be given functions of t. Denote 

cos a sin a cos (3 sin a sin j3 
T—\ — sina cos a cos (3 cosa sin/? | (3.1) 
— sin /? cos 13 



and 



Then 



a' (3' sin a 

Q = ( -a' (3'cosa \ . (3.2) 

—(3' sin a —/3' cos a 

cos a — sin a 

rp-i^rpt^ j sin a cos/? cos a cos/? -sin/? | (3.3) 
sin a sin /? cos a sin /? cos /? 



and 



|(T) = gr. (3.4) 



Define the moving frames: 







/ X 


y 






2 I 







U \ ( u 

V \=t[ V \ , \ y \ =T\ y \ . (3.5) 
I \w 

Note 

A = a^ + a^ + a2 = 5| + a^ + 9|, ud^ + vdy + zd,^udx + vdy + wdz, (3.6) 

«. + ^, + ^i;.=W;, + Vy + >V2, 1^ j;* j =g 1^ 3; j , (3.7) 

dx\ ( d,\ I dt{U) \ ( U\ ( ut\ 

dy \=T\ dy \, dt{V) \=Q\ V \+T\ V, . (3.8) 
dz I \d, J V J \W J \wt J 

Write W, V, W,p as functions in t, X, y, Z. Set 

Ri^Ut + a'{yUx - XUy - V) + P'{ZUx - XUz - W) sin a 

+P'{ZUy - yUz) cos a + + VUy + - J/A(W) , (3.9) 

i?2 = Vt + a'(3^VAf - XVy +U)+ (3'{ZVx - XVz) sin a 

+(3'{ZVy - yVz - W) cos a + UVx + VVy + WV^ - j/A(V), (3.10) 

13 



i?3 = + a'{yWx - XWy) + l3'{ZWx - XWz + U) sin a 

+P'(ZWy - yWz + V) cos a + UWx + VWy + WWz - i^A(W) , (3.11) 

Then Navier-Stokes equations (1.1)-(1.4) become 

Ri + -p^^O, R2 + -Py^0, Rs + -p^^O, (3.12) 
P P ^ P 

Ux + Vy + Wz = 0. (3.13) 

Instead of solving the equations in (3.12), we will first solve the following compatibility 
equations: 

dyiRi) = d;,{R2), dz{Ri) = d;,{Rs), dz{R2) = dy{R^) (3.14) 

for V, W, and then find p from the equations in (3.12). 

Let / be a function in t,y, Z such that dy{f) — d%{f) — 0, and let 0, be functions 
in t, X. Suppose that 7 is a function of t. Assume 

U^f-2iX, V^(l) + iy, W^i/J + ^'Z. (3.15) 

Then 

Ri = ft- 2/ A" - a'{3-f'y + Xfy + 0) - P'{3iZ + Xfz + i/j) sin a 

+(3\Zfy - yfz) cosa - 27'(/ - 27' A") + fy{^ + 7'^;) + /^(V^ + 7'Z), (3.16) 

= 0t + i'y + a'{y4)x - 37' A" + /) + l3'Z4)x sin a - jS'il^ cos a 

+(/ - 2-i'X)<i)x + 70 + {ify - ^te, (3-17) 

Rz = ^t + l"Z + a'yxl^x + ^'{Zil^x-^iX + f)sina-v'il:xx 

+/?Vcosa + (/ - 27'A')V;;t + 7'(V' + V^^)- (3.18) 

Now (3.14) becomes 

(t)tx + {a'y + sin a + f )(l)xx - P'i'x cos a - 2j'{X4)x)x 

+i(t)x - v4>xxx = fty - P'fz cos a - j'fy, (3.19) 

i^tx + (a'y + P'Z sin a + f)ipxx - ^i^xxx + P'(t>x cos a 

-2{'y'XiJx)x + ii^x = ftz + P'fy cos a - ifz, (3.20) 
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a'fz + ((3' sin a + fz)(t>x = («' + fy)i^x + fy sin a. (3.21) 
By (3.19) and (3.20), we take 

f ^ -a'y - l3'Zs\na (3.22) 

modulo the transformations of type Ti^ in (1.8) and (1.9). Note that (3.21) is imphed by 
(3.22). Integrating (3.19) and (3.20), we obtain 

0t — 2'^' X(j)x + 7V — ^(tixx — /3'ipcosa = [(/3')^sina cos a + a'^' — a']X + (3.23) 

tfjt — 2j'Xipx + tV ~ ^i^xx + cos a 
= -[(/?'sinQ;)' + Q;'/3'cosQ;-7'/3'sina]A' + /32, (3.24) 

where /3i and are arbitrary functions of t. To solve the above problem, we write 

/3'=^, 7=iln//' (3.25) 

cos a 4 

and set 

y ^ y sm(^ cost/? J \ W J 



dt. (3.27) 



71 A /"^(^ ^"''^ '^^'^ ^ ( (^')'tana + ^ 

72 / J W\ -sinv? cos(/J y _((^'tana)' - aV' + ^ 

Then (3.23) and (3.24) are equivalent to: 

a" 

0t - ^,^(t>x - i^(t>xx = ii^'X + 99;, (3.28) 
u" 

ipt - - y^^xx = X + (^^, (3.29) 

where and (/72 are arbitrary functions of t. Thus 

m 

= 7iA/i^'^ + ^1 + E cZi,ie'^(«t'-''M)M+«Mvy^ sin(6i,,(2ai,,// + ^i^/A') + q,i), (3.30) 

1=1 

^ = 72 + <^2 + 5^ cii,,e'^("t-fU'^+«2,.V^-^sin(62,.(2a2,.// + ^ X) + cs,.), (3.31) 

s=l 

where ai^j, 02,5, 62,3) Ci,i, C2,s, and 0^2,5 are real constants. According (3.26), we have 
= y;;7(^^cos(/p + 72sin^)A' + ai + ^^c^i,,e'^K-''l>+«Mx/^^ 

n 

X sin(6i,,(2ai,, + ^X) + q,0 + ^ ^ rf,^^e^(«i.-''f.J/^+«.,.V^-^ 

v/^ s=l 

X sin(62,.(2a2,. + ^f]j X) + C2,.), (3.32) 
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n 

^ ^ s=l 

X sm(62,.(2a2,. + y^X) + ca,,), (3.33) 



where ai and (T2 are arbitrary functions of t. 
To find the pressure p, we recalculate 

= {{ip'yy-2ip''il;)tana-2a'(l)-a"y-{(p" + a'(p'{l + sec'^a))Z 

,"iay + ^'Z tanc) ^ ^^^,^2 ^ ^ , (3.34) 



„ (4/'-3(/)2)j; , ,,2^, , , , (a'li" + a"ii')X 

^3 - ^tana + a, 

-{ip'yztan^a-{ip" + a'ip'{l + sec'^a))X. (3.36) 

Thus 

p = p{(a"^ + (/tana + «y(l + sec^«))Z);^+^"^"'^ + ^'^'""")'^ 



+ 2 J{a'(j)+(p'i)tana)dA:} (3.37) 



2 

modulo the transformation in (1.10). In summary, we have: 

Theorem 3.1. Let a,ip, fj,,ai,a2 be functions oft with /i' > 0. Take real constants 
{0'i,i, ci2,s, h,i, b2,s, ci,i, C2,s, di^i, d2,s M = 1, m; s = 1, n}. Denote 

/3=/^ (3.38) 
J COS a 

and define 71,72 by (3.27). Take the notations X,y,Z given in (3.1) and (3.5). In 
terms of the functions (p in (3.32) and ip in (3.33), we have the following solution of the 
Navier-Stokes equations (1.1)-(1.4): 



u 



— — -)_ Q,'y _|_ (p'Ztajia] cos a — ( \ sin a, (3.39) 

\2ii' J \ 4/1' J 
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V — 



oly — Lp' Z tan a ) sin a cos (3 



+ (^0+^^ cosa cos/3- + ^ ) sin/3, (3.40) 



w — 



a'y — if'Z tan o; ) sin a sin P 



+ (^0 + ^^ cosa sin/3+ + cos/? (3.41) 

and p is given in (3.37). 

Remark 3.2. We can use Fourier expansion to solve the system (3.28) and (3.29) 
for 0(yU, \/Ji!X) and ^(/U, \/Ji! X) with given 0(0, ^/Ji' X) and ^(0, ^/Ji!X). In this way, we 
can obtain discontinuous solutions of Navier-Stokes equations (1.1)- (1.4), which may be 
useful in studying shock waves. 

Let /, g, h be functions of t, X , y, Z that are linear in X ,y, Z and fA:+gy + hz = 0. 
Based on our experience in [X2], we assume 



U ^ f -6iyX-\ V^g-6iyyX-\ W ^ h. (3.42) 



Then 



Ri = ft + fh + fy9 + fzh-6uf;,X-' + a'{yf;,-Xfy-g) 
+f3'{Zfx - Xfz - h) sin a + f3'{Zfy - yfz) cos a 

+6u{f - yfy + 2a'y + P'Z sin a)X-^ - 2Au'^X-\ (3.43) 



R2 = gt + f9x + 99y + 9zh + (y-'{y9x - ^9y + /) + P'{Zg^ - Xg^) sin a 
-Qv{a' + g^)X-^ - 6u{g + (3' Z cos a-a!X + y9y)X~^ 

+P'(Zgy - g^y - h) cos a + I2vy{f + a'y + {3' Z sin a) (3.44) 

R'i = ht + fhx + ghy + hhz + a'iyhx - Xhy) + [3'{Zhx - Xhz + /) sin a 

+P'{Zhy - yhz + g) cos a - Qu{h^ + (5' sin a)X-^ - 6u{hy + (3' cos a)yX^'^. (3.45) 

According the negative powers of X in (3.14), we take 

f ^r^X -a'y - 13'Zsma, g ^ a'X + -fy - P' Z cosa, (3.46) 

h= -{P'Xsma + p'y cos a + 2-fZ) (3.47) 
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modulo the transformations of type in (1.8) and (1.9) for some function 7 of t. With the 
above data, we have: 

+(3(/?')^ sin a cos q; - a" - 2a' + (4/?'7 - /?")Z sin a, (3.48) 

i?2 = {i + -{a'f + ^{(3'fcos^a)y -l2pa'X-^ 

+(a" + 2a'-i + sin a cos a) A" + (4/?'7 - I3")Z cos a, (3.49) 

i?3 = (472 - 27' - {l3'f)Z + (4/3'7 - /3") (A" sin a + cos a). (3.50) 
By (3.48)- (3.50), (3.14) is now equivalent to 



-Q;"-2a'7 = Q;" + 2Q;'7=^7 = -;^. (3.51) 

2a 



Thus 



-—X-a'y-P'Zsma-QvX"\ V = a'X-^y-fi'Z cosa-6uyX-\ (3.52) 



a" 



W= — Z-/3'A'sina-/3'3^cosa (3.53) 

by (3.42), (3.46), (3.47) and (3.51). Moreover, (3.12) and (3.48)-(3.50) imply 

^ {2a'a'" + 4{a')^ - 3{a" f) {X^ + y^) 3{f5')\X^ sin^ a + y^ cos^ a) 
P - 8(aO' 2 

+ 12i/(a'3^A'"i - uX-^) + (/5" - 4/5'7)Z(A' sin a + 3^ cos a) 

./o/N2^., . (Q;W" + (a')'(/3')'-2(a'7)^% /or.x 

-3(/3' 'YJV sin a cos a + ^ ^ Z L ^ — — — }• (3-54) 

2[a')'^ 

By (3.3) and (3.5), we have the following theorem: 

Theorem 3.3. Let a and P be functions of t with a' 7^ 0. In terms of the notations 
X,y,Z given in (3.1) and (3.5), we have the following solution of Navier- Stokes equations 
(1.1)-(1.4): 

u=\ — - + 6i'yX ^ J (3^ sin q; — A* COS a) — a' (A" sin a + 3^ COS a), (3.55) 
\2a' J 

f a" _\ 

V = — [ — - + QuyX ^ ) (A sin a + 3^ cos a) cos (3 + a'(X cos a — y sin a) cos P 
\2a' J 

-P'Z cos P + (^P'X sin a + p'y cos a - ^Z^ sin /3, (3.56) 
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w — —\- VQi^yX \ iX sina + y cos a) sm 13 + a\X cosa — y sm a) sin 13 

\2a' J 

-j3'Z sin 13 + {^Z - 13' X sin a - f3'y cos cos ^ (3.57) 

andp is given in (3.54)- The above solution blows up at any point on the following rotating 
plane with a line deleted: 

{{x, y,z) eM.'^ \ X cos q; + y sin q; cos /3 + zsina sin P — 0, 

—X sin a + y cos a cos (3 + z cos a sin /3 7^ 0} . (3.58) 



We remark that the above solution may be apphed to study turbulence. 

4 Moving-Frame Approach II 

Motivated from the solution ip in (2.27) of the equation (2.24), we will solve Navier-Stokes 
equations by ansatzcs with given irrational functions under the moving frames in (3.5). 
First we rewrite (3.9)-(3.11): 

Ri=Ut + {a'y + 13' Z sin a + U)Ux + (V - a'X + (3' Z cos a)Uy 

+ (>V - 13' {Xsma + y cos a))Uz - a'V - 13' W sin a - i^A{U), (4.1) 

i?2 = Vt + (a'y + 13' Z sin a + U)Vx + {V - a'X + 13' Z cos a)Vy 

+(>V - 13' {X sin a + 3^ cos a))Vz + a'U - p'W cos a - i^A{V) , (4.2) 

Rs^yVt + {a'y + 13' Z sin a + U)Wx + {V - a'X + 13' Z cos a)Wy 

+(W - 13' {X sin q; + cos a))Wz + 13' {U sin a + V cos a) - vA{W) . (4.3) 

Let cci, 7 be functions of and let a, h be real numbers. Set 

^^ = e^ry+PiZ _^^-c^^y-PrZ^ ^o = e"i^+'3i^ + oe-"i^-^^^, (4.4) 

a = sin(ai3; + l3iZ), (i = cos(ai3^ + ^Z), (4.5) 

00 = e-y^ - be-^^ , Co = e^-^ + be'^^ , (4.6) 

(Pi = sm{^X), ^i = cos(7A'), Ai = 4 + a|. (4.7) 
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Suppose that / and h are functions in t,y,Z. According to (3.46) and (3.47), we 
assume 

U - -a'y - ^'Zsina- {fy + hz)X - {a,a + AT)Cr</'., (4.8) 

V = a'X - P'Zcosa + f + a-f^rips, W = sina + y cos a) + h + T'j^ri/Js- (4.9) 
By (4.1)-(4.3), we have 

R, = -{a^a + /?ir)'Cr0. - {^la + Pir)[{-ina[y + f3[Z)^r<P, + 

-ifyt + hzt)X + {{fy + hz)X + (aia + /3ir)Cr.0s)(/y + hz + 7(aia + /3ir)Cr^.) 

—a'{f + olX — (3'Z cos a + ^cr^rfps) — P'iP'i^ sin a + 3^ cos a) + h + ^r^rips) sin a 

-(/ + 7c^6V's)(«' + (/yy + hyz)X + {-lyaiiaia + Pit)CAs) -{h + IT^M 

x(/3'sina + {fyz + hzz)X + (-l)'-/5i(aia + Pir)^r(ps) + ^/{Ai(/y + /i^)^' 

+{aia + PiT)[{-iy{al + Z^^) + (-i)^f ]C,c/*,} - q;"3^ - (/3'sina)'Z 

= {(7(/y + - l')XCrA - i-ina[y + (3[Z + a,f + /3i/i)e.0j(Q;i(T + /3ir) 

+{(«ia + /3it)[(-1)V72 + (-l)V(a? + Z^^) + + hz] - {a^a + r/3i)'}C.</>. 

-7{2(aa' + r/3' sin a) + [(r{fyy + /iy^) + r{fyz + Z^^^)] ^'j^.V. - (/yt + hzt)X 

+{fy + hzfX - f{a' + ifyy + hyz)X) - sin a + {fyz + hzz)X) 

-a\f + olX - P'Z cos a) - p'{P'{X sin a + 3^ cos a) + h) sin a - a"y 

-iP' sin ayz + uAii fy + hz)X + 7(ai(T + i3^Tf{5r,i + 4a5,,o)</'s^s, (4.10) 

i?2 = c^''^" - (/3'cos q;)'Z + /t + (7<j)'erV'. + l^{{<^'iy + /3i2:)C.V'. + (-l)^7''^er<^.) 
-u'lay + sin a + (/y + hz)X + (ai(T + /3iT)Cr0s] - I3'[I3\X suia + y cos a) 

+/i + 7rerV',] cos a - [(/y + /i2)A' + (aid + /5ir)Cr0s](a' + (-l)'7'cr60.) 
+(/ + l(^^r'tps){fy + aaaCrtps) + {h + 'yr^rips){fz - P' cos a + Pi^yaCri^s) 
-^[Ai(/) + ia{{-inal + /3,^) + (-l)^f )]eV. 

- + + 7a(«;3^ + P[Z + ai/ + A/i)Cr-^s + (-1)^7^(7 - 7(/y + hz))Xir^, 
+H^fy + ^/^ - Mi-^Yiotl + /32) + (-1)^72] - 2r/?'cosa] + (7^)'}^^^. 
-2a'{aia + /3iT)Cr<^. - (/3'cosq;)'Z - a'fa':^ + sin a + 2(/y + hz)X] 
-P'[P'{X sin q; + cos a) + h\ cos a + //y + h{fz - P' cos a) - uAi{f) 
+7V(aia + /3ir) (46^3,0 + 4,i)CrCr, (4.11) 

i?3 = iP' sin a)'A' + (/?' cos a)y + ht + (r7)'e.Vs + r^[(a[y + P[Z)C^, 

+ {-irj'XCrC^s] - [{fy + hz)X + («ia + P,T)Cr(p.]{P'sma + {-l)W^rC^s) 
+ {f + (^'y^ri>s){P' COS a + hy + aiT'y(r'ips) + {h + T'y^ri>s){hz + PlT^CAs) 
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-P'ia'y + P'Zsma + (fy + hz)X + {a^a + (5iT)Cr(t>s) sina + /?'(«' A' - (3' Z cos a 
+f + c^l^ri^s) cos« - iy[Ai{h) + 7r((-l)^'(a? + Pi) + (-1)^7')]?^^. 
= ^T{a[y + P[Z + a,f + mCr^s + {(r7)' - ^^7^ + + (-1)^7'] 
+7(2/3'crcosQ; + cr/iy + r/i2)Kr'0s - 2/3'(Q;ia- + /3iT)Cr0s sina + (-1)*7t(7' 
-lify + /i2))'^er0s + (/?' sin a)'X + (/?' cos a)'y + ht - (3'{fy + /i^) A" sin a 
+f{P' cos a + hy) + /i/i^ — l3'{a'y + sin a + {fy + hz)X) sin a + /3'(q;'A' 

-/3'Zcosa + /) cos a - uAi{h) + -f'^T{aia + ^iT){4bSs,o + Ss,i)^rCv (4.12) 

By the coefficients of ^ri^s in the equation dy{Ri) = dx{R2), we have 

7V = {-lY+'+'a,{aia + [a{fyy + hyz) + T{fyz + /l22)]3; = 0. (4.13) 

Moreover, the coefficients of C^^s in the equation dy{Ri) — dx{R2) suggest 

{fy + hz)y = 0, (4.14) 

which imphes the second equation in (4.13). According the coefficients of ^r4>s in the 
equation dy{Ri) = dx{R2), we get 

aPihy = ai{Tfz - 2t 13' cos a). (4-15) 

Furthermore, the coefficients of C^V^s in the equation dy{Ri) — dx{R2) yield 

ai(3' svna = a' Pi. (4.16) 

Symmetrically, we have (4.16), 

7V = (-1)'-+-+Vi(«i(^ + At), {fy + hz)z = (4.17) 

and 

fz^hy + 2(5' cos a. (4.18) 
By the first equation in (4.13) and (4.17), we have 

cr/3i = Tai. (4.19) 

Then (4.15) is implied by (4.18) and (4.19). Note that the equations of the coefficients 
iri^s, Cr'fps-, Cr(l>s and Cr-0s in dz{R2) — dy{Rs) are implied by (4.16), (4.18) and (4.19). 
According to (4.14) and the second equation in (4.17), 

fy + hz^li, (4.20) 
21 



a function of t. Under the conditions in (4.16), the first equation in (4.17), and (4.18)- 
(4.20), dy{Ri) = dxiR2) becomes 

a'hz — P'hy sin a — a'' , (4-21) 

dz{Ri) — dx{R3) is equivalent to 

(3'hz sin a + a'hy = (3'^i sin a — {(5' sin a)' — 2a' [3' cos a (4.22) 

and dz{R2) = dy{Rs) says 

iffy + hfz)z = ifhy + hhz)y + 2/3'7i cos a. (4.23) 

By (4.18) and (4.20)-(4.22), we assume fy, fz, hy and hz are functions of t. Then (4.23) 
can be written as 

ify + hz)fz = ify + hz)hy + 2f3'j, cos a, (4.24) 
which is imphed by (4.18) and (4.20). Solving (4.21) and (4.22), we get 



a'P'^i sin a — {a'/3' sin a)' — 2(q:')^/3' cos a /a nr.\ 



a'a" + {p'^ji sin^ a - {(3' sin a) {(3' sin a)' - a\(3'f sin 2a 
{a'Y + {(3'fsm^a 

Moreover, 

7i(a')^ - a'a" + (/3^sina)(/?^sina) ^ + a'(/?')^ sin 2a 
(Q;')^ + (/?')^sin2a 

by (4.20) and (4.26), and 



^2 = /.A9 , /^A2.^„2,_ • (4-26) 



a;'/3'7i sin a — {a' (3' sin a)' + 2(/3')^ sin^ a cos a 
(Q;')^ + (/?')^sin2a 



u; 13111 oi — \^u; fj 13111 u; ^ -p ^(^/O/ ) isiii u; i^uis oi / ,1 \ 



by (4.18) and (4.25). With above data, we take 

f = fyy + fzZ, h = hyy + hzZ (4.29) 

by the transformations of the type in (1.8) and (1.9). Furthermore, (4.17), (4.19) and the 
first equation in (4.17) yield r + s + 1 e 2Z, 



ax = <^a', 7 = ±^-y(a')2 + (/3')^sin2«, (4-30) 
(3\ — (pf3' sin a, a — ^a' , r = ///J'sino;. (4-31) 
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In particular, a,P,^i,(p and are arbitrary functions of t. According (4.9)-(4.11), the 
pressure 



p = p{7/x^-^[(y - jji)XCr(l)s - ii^ayy + (^/3'sin ayz + ifia'f + /3'/isina))e.V.] 
+(-1) V-^[(7//)' - ^iicp'^-'Kr^s + 2//((«')' + iP'Y sin^ c^%^s + 2a' fX 

(q:')- + sill- n ; -7? , , r. , ^, , , n 

+ — ^ ^ ^ ^^"2 + 2(3' hX sin a + [(/?' sin a)' - a'/?' cos a]XZ 

+ sin 2« + a"^ A-Jt^ - l7VV"'[(^ni + 4a5.,o)</'^ + (4W.,o + 

+[(/?' cos a)' + a'/3'sina - /^^ - /^/^ - hyh^\yZ + W - /I - 4 ^2 
^ K)^ + (/30^cosa-/y,-/^-/.^ ^,,^ 

modulo the transformation in (1.10). By (3.3) and (3.5), we have the following theorem: 

Theorem 4.1. Let a,/?, 71, <^ and fj, be arbitrary functions oft such that 7^ and 
{a')^ + (/3')^sin^a ^ 0. Take any real constants a and b. The notations X,y and Z are 
defined in (3.5) via (3.1), and the notations $,r,Cr,(t>r and if^r are defined in (4.4)-(4-7) 
with ai, (3i and 7 given in (4-30) and (4-31). Moreover, fy, fz,hy, hz and f, h are given 
in (4.25)-(4.29). Assume (r, s) G {(0, 1), (1, 0)}. We have the following solution of the 
Navier-Stokes equations (1.1)-(1.4): 

u — — a' {X iiuia + y cos a) — {f + iia''^^r^s) ^^"^ Oi 

-{-^^X + ipii{{a'f + {P'f srn^ a)Cr(t>s) cosa, (4.33) 



V — {f cos a — P'Z) cos /? — {a' sin a cos P + (3' cos a sin P)y 

— {jiX + ipn{{a'Y + (/3')^sin^ a)Cr4's) sin a cos/3 — hsmf3 

+ (Q;'cosa; cos /3 — sin sin (3){X + j/i^r'ips), (4.34) 

w — {(3' cos a cos (3 — a' sin a sin Pi)y + (/ cos a — (3' Z) sin (3 

— (71 A' + (pii{{a'Y + (/3')^ sin^ a)C,r4>s) sin a sin /? + hcos^ 

+ (a' cos a sin /3 + sin a cos/3)(A' + 7yu^r'0s) (4.35) 

and p is given in (4-32). The above solution is globally analytic in x,y,z. 

Let 7i, 72 be functions of t and let a, b, c be real numbers. Denote 

00 = eTi^+T2^ - 01 = sin(7i:y^ + 72Z), (4.36) 
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^0 = ei^y+'^^z ^ ^^-71^-722^ ^ cos(7i3^ + , (4.37) 

^0 = 6e^i^+^^^ - ce-^i^-^^^, a = csin(7i3; + 72Z + 6), (4.38) 
(0 = &e^iy+722 + cg-7i3^-72^ ^ ^ c cos(7i3^ + 72Z + 6) . (4.39) 

Suppose that a, r are functions of t and /, k, h are functions in t, X, y, Z such that h and 
5^ are hnear in X,Y,Z and 

+ A;^; + /i^ = 0. (4.40) 

Motivated from the above solution, we consider the solution of the form: 

U = -a'y - p'Zsma + f - {jI + j2)HrX + ai/jrX^), (4.41) 

V^a'X- (3'Zcosa + k + 71 (r^^ + 2a(j)rX), (4.42) 
Vy = ( A" sin a + cos a) + /i + 72 (r^^ + 2f70^ A") . (4.43) 
For convenience of computation, we denote 

7 = 71+72, r^f-f.X Ai = 4 + a|. (4.44) 

Now (4.1) becomes 

R, = -a"y - iP' sin a)'Z + /, - (-l)^7(7i3^ + 72^^) (r^r^ + a<t>rX^) 
+((-l)V7V - M')CrX + (/ - jirCrX + aiPrX'mfx - lirCr + 2<J'll^rX)) 

+ {k + 71 «r + '2a4>rX))[fy - 2a' - (-l)^77i(rC,A + a^.A^)] - z/Ai(/) 

+ {h + 72(Ter + 2(J0,A))[/2 - 2(5' Sma- {-iy-i-i2{TirX + (70,^2)] + 2i/7(7^, 

-a'(a'A'-/?'ZcosQ;) - (/?')^(A'sinQ; + ^^cosa) sina + ((-1)^7^ - (7(7)')^'^'^^ 
- -((a)^ + (/3')^sin^a)A - (a" + 2-^(/5')^ sin2a)3^ + (a'/3'cosQ; - (/3'sina)')2^ 
+7^[r^(465o,r + c5i,r)cA' + 3ar(25o,r(a& + c) + (51,^0 008 6)^^ + 2cT^(4a(5o,,. + (5i,r)A'^] 
-(-l)''7(7i3^ + 72^ + ^71 + ^72)(re.A + cr^.A^) + //^ + kfy + 
+((-l)V7V - (7(7)' - 37a/^)V;.A'2 + ^2^^^'^ - Ai(/)) - 7x7* C. 
-[((7r)' + 27r/;, - (-1)V7V)C. + 27arV,]A' + ft 

+ (7i(/y - 2a') + 72(/2 - 2/3'sina))(r^, + 2a0,A). (4.45) 
To solve (3.14), we assume 

i^y + 7^Z + A;7i + /i72 = (4.46) 

and 

(_1)V7V - (7a)' - 37a/;, = 0, (4.47) 
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Moreover, (4.2) and (4.3) become 

R2 = a"X - iP'cosayZ + ((7ir)' - (-l)V77ir)e. + 2((7ia)' - (-l)V77ia)0,A' 
+kt + {-f[y + 72^^)71 (-rCr + 2a7PrX) + (/ - lirCrX + a^rX^)){2a' + kx + 271 (70,) 
+ {k + 7i(r^, + 2a(prX)){ky + 7i(rCr + 2(7?/;, A')) - {l3'f{Xsma + ycosa) cos a 
-a{a'y + /3'Zsina) + {h + 72 (r^, + 2(7(/),A'))(A;2 - 2(3' cos a + 7172 + 2a?/', A*)) 
= {a" - 2-\p'f sin 2q; + fx{2a' + kx))X - {{a')^ + {(5'f cos^ a)y + kt + kky 
+[^(71^ + 72(^2 - 2/3' cos ck)) + (7ir)' - (-1)V77it]^, - {{13' cos a)' + q;'/3' sin q;)Z 
+7ct(27ict(/), - 2q' - kx)A'^'^ + f*{2a' + kx + 271 (70,) + h{kz - 2/3' cosct) 
+77i^'^rCr + {277i(7r^,V^, + 2[(7i(7)' - (7(7i(/i^ + (-1)^7) 

+72(2/3' cos a - kz))\(t>r - 1t{2oI + A;;t)Cr}'^ (4.48) 

i?3 = (/3' sin a)'X + (/3' cos c^)':^ + (72t)'6 + 2(72(7)'0,A' - (-1) V7|(Te, + (70,A') 
+ (7l>' + 722)72(rCr + 2c7?/;,A') + (/ - 7(rC,A' + c7V',A'2))(2/3' sin a + /i;, + 272C^0r) 
+(A; + 7i(t^, + 2(70,A'))(2/3'cosct + hy + 7i72(rC, + 2(7^/;,^*) - {i^'fZ + /ij 
+q;'/3'(A'cosck - J^sina) + (/i + 72 (r^, + 2(70,A'))(/i2 + 7^(tC, + 2(7-0, A")) 
= [(/3' sin ck)' + ct'/3' cos ct + /;t(2/3' sin a + ^a^)] + [(/3' cos a)' - q;'/3' sin a\y 

+ [(72^)' + (71(2/3' cos a + hy) + 72/1^ - (-l)V772)r]e, + {2^^2T<jir^r + 2[(72C7)' 

-72(7(/cy + (-1)V7) + 71(7(2/3' COS a + /iy)]0, - 7t(2/3' sincK + /i^)(^,}A' 
+/*(2/3' sin q; + + 272(70,) + k{2^' cos a + hy) + ht + hhz + 772T^CrCr 
+7(77/^,(272(70, - 2/3' sin a - hx)X^ - {(3'fZ. (4.49) 

By the coefficients of A"^ in dz{R2) — dy{R3), we have: 



According to (4.46), 

kxli + hx^2 = 0, 7i + ^iky + 72/z.y = 0, 72 + ^^kz + 72^2 = 0. (4.51) 
Solving (4.50) and the first equation in (4.51), we obtain 



72(2a' + kx) = 7i(2/3'sina + hx). 



(4.50) 



kx = 27 S2(/3'7i sin a - q;'72), 



hx = -27 ^7i(/3'7isinQ; - q;'72). 



(4.52) 



Moreover, the coefficients of X in dz{R2) — dy{Ri) give 



7^72 - 7i72 + 7i72(^y - hz) + 72^2: - llhy - 27/3' cos a = 0. 



(4.53) 
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By (4.51), the above equation can be rewritten as 

kz-hy^2p'cosa. (4.54) 

Furthermore, (4.50) and the coefficients of in dz{R2) = dy{R'^) show that / is a 
function of t and 7i3^ + 72^^. According to the coefficients of X in dy{Ri) = dx{R2) and 
dziRi) = dxiRs), we take 

/* = (fdr + aw(f)r + «!, (4.55) 
where (p and cui are functions of t, and 

^ = 7i3^ + 72Z, ^0 ^ bie^ - cie-^ , t?i = Ci sin(ro + 61) (4.56) 

for 61, ci e ]R. 

Now the coefficients oi X in dy{Ri) = dx{R2) and dx{Ri) = d^iRs) give 

r[2ai7(7V'. + iM' + 2^rfx - (-1)^7^)^] = 0, T^dy,dz. (4.57) 

Moreover, the coefficients of X^ in dy{Ri) = dx{R2) and dx{Ri) = dx{Rz) yield 

[(/;, - {-iriv)v + v')'&r + ((/^ - (-l)V7)a + a')^0. - Q;i7<.]y 
= 2[(7iCT)' - a(7i(/i2 + (-1)V7) - 72/^^)]^^ 

+2a" + 2a' fx + A;;^,^ + - kxhz, (4.58) 

[(/;, - {-iriv)v + (^0^. + ((/;t - (-l)V7)a + a')^0. - a^irQz 
= 2(/?' sin a)' + 2(3' fx sin a + /i^-t - ^/^a- + kxkz 

+2[(72(7)' - 72Ct(% + (-1)V7) + 7i(xfc2](/>.. (4.59) 

Thus we have: 

2a" + 2a' fx + A:;(^t + - kxhz = 0, (4.60) 

2{I3' sin a)' + 2/?'/^^ sin a + hxt- kyhx + kxkz = 0. (4.61) 
For simphcity, we only consider two special cases a follows. 

Case 1. 'dr = Cr, = 0, 71 = a'/j, and 72 = f5' iJ,sma. 

In this case, 

kx = hx = (4.62) 
by (4.52). According to (4.60) and (4.61), we have 

a" 

p'sma = da', fx = • (4.63) 

a' 
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Moreover, (4.51) becomes 

ky + dhy = -i^, kz + dhz = -d^-^. 

According to (4.40) and (4.54), 

a" 

hz = ky, hy = kz — 2(3' cos a. 

a' 

Substituting (4.65) into (4.64), we obtain 

ky + dkz = 2dP'cOSa- kz - dky = _^(^f^Z±^ 

fj,a' fj,a' 
For convenience of computation, we write 



From (4.66), 



' 2d^a' cot a - . 



By (4.65), 



1 + d^ \ (3[ a' 

d /„ , 1%' a"\ 



hy = 

Furthermore, (4.57) becomes 

(7r)' + 2-frf;, - (-1)V7V = ^ 7r = (aO'e^"'^""^'^"'^'^'- 

So 

(a;')^e("^)''''(^+'^')'^i 



il + d^)/3[ ■ 
Note that (4.58) and (4.59) are imphed by 

,„ , . , „ aV-{a"+{-iru{l + d')a'P[)^ 
ifx-i-l) + ^ -ai7r = 0^ai = (^;)3e(-i)-.(i+<i^)/3i " 

Observe 

^ = -^A'-a(j; + dZ) + (^.)3,(-i)^.(i+.^)/.. 

+ (^_(«')2e(-i)Mi+<i^)/3i;t')C., 
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^ (2(fa'{y + dZ) cota + {{d^ - l)y - 2dZ)^ + {d^y - dZ)^] 
\- d"^ \ ipi a' J 



1 + ^2 

+a'{X - dZ cot a) + ^ \^ ^ — ir 



(l + d2)v^/3i 



(4.76) 



1 / R' ' 

w = ,^^U{i-d^)z-2dyy-^^{z-dyi^-2d{y + dz)olcot 



a 



^dol{X + :ycota) + ^ ^ 



(4.77) 



Moreover, 

i?i = (465o,r + c5i,,)c(a')'e(-i)^'^(i+'^')^^ [(a')'e(-i)^^(^+'^')^^ A' - ^] 



(2(a")2 - a'")A' 2(a')^e 



,/\3.(-l)V(l+d2)/3i 



+o?((a')^ cot OL - a")Z - {a" + {da'f cot a):^, 



- (a') (1 + 



(4.78) 



R2 = 7i 



(T+n^^ 



-(a" + (da')^ cot a)X + [A;^* - {a'f{l + rf^ csc^ a)]:V^ + -(A;^ + /i^)^ 
+\{hz - /3'cos«)' - d(a')']^ + 2«7* - 2(Q;')'e(~')''"^'+'^')^iCrA', 



(4.79) 



^3 = 72 



(l + d2)(/3i)2 



+d((a')^ CSC q; - cl')X + ((/ly + cos a)' - d(a')^) + ^(^^^ + /i^)^ 
Hhzt - {da' f esc'' a)Z + 2da7* - 2d{a'fe^-'-'^^''^^^'^"^^'CrX. 
By (3.12), 

p = p{^^ ^= ^ + ^^(y + d' csc^ a(3;2 + z') + 2rf3;Z) 

ji\2 ^ti'\ v2 ^y^,/'' ^'„/'^2 



(4.80) 



^ {2{a"Y - a"')X' ^ a! a!" - {ot 



+ 



) 222 (Q;')''e(-^)'"2''(^+'^ ^'^i,^^ 

2(a')' ' 2(l + d2)«-^2rf3^^"^^ 2(1 + d2)/3i ~ 

l)^a'(a'/9r - 2a"/?0Cr 2a'[(a" + (-l)V(l + d')a' ii'^)^ - aV'lC^ + 



(l + (i2)(/3;)2e(-i)'-+^-(i+<i^)/3i 
1 



+ 



[ 2d'cl{y ^dZcoXa) + 



(Q;')3e(-i)''''(i+'^'')^i 
^'^{{^ - l)y - 2dZ) a"{d'y - dZ) \ 



+ 



2(l + d2)2 

(5'^{{1- d')Z -2dy) a"{Z-dy) 



2f3[ 



+ 



a' 



2(5[ 



+ 



a' 



2d{y + dZ)a' cot a 



[a'f{l + d'')X'' 



d 



[a" cot a - {a'Y esc' a)((l - rf'););^ + 2diZ' - y')) + 



' 1 + d^' V / . V - // . 4(l + rf2)(/3;)2 

X (4d:y^Z + (1 - d2)(:y^2 - Z')) + (a" + {da'f cot a)A'3^ + (i(a" - {a'f cot a)XZ 
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(4.81) 



modulo the transformation in (1.10). 

By (3.3) and (3.5), we have the following theorem: 

Theorem 4.2. Let a,Pi,(p be functions oft and let b,c,d be real constants. Denote 

/? = (iln I csca — cotcul (4.82) 

(so the first equation in (4-63) holds). Define the moving frame X, y and Z by (3.1) and 
(3.5), and 

^ be^i{y+dz) _ ce-^^(y+'i^\ ^^^csm[^/p'i{y + dZ) + b], (4.83) 

Co = 6e^(^+'^^) + ce-^(^+'^^) , Ci = c cos[ V^C^ + dZ) + b] . (4.84) 
For r — 0,1, we have the following solution of the Navier-Stokes equations (1.1)-(1.4).' 



U 



a'i^' - {a' + (-l)V(l + d')a'(3[)i^ _ , , o^.^y.^^+a^^^ . _ 
-^-^ (^d-'a'iy + dZ) cota + ((c/2 - l)y - 2dZ)^ + {d'y - dZ)^^ si 



cos a 



+a'{X sin o; — (^ + 2dZ) cos a) H ;-; — — 7= — sin a, 



sm a 



(4.85) 



(Q;')3e(-i)''Ki+rf2)/3i ^^'^ 



sin q; cos P 



^2da'{dcosa cos/? — sin/3)(3^ + dZ) cota ^ Q;''((icosQ; cos/3 + sin/?)(d;V — Z) 



1 + d^ 



(1 + d^)a' 



[((d^ - l)y - 2dZ) cos a cos /? - ((1 - d'^)Z - 2dy) smp]p[' _ a" 

2{l + d^)(3[ a' ^ 

+q;'A'(cosq; cos/3 — (isin/3) — a'3^(cos a cos /? — cos a sin/3) 



, ,^ ^ (a')^e(-^)'^(^+'^')'3i(cosa cos/3-dsin/3), 
-da Z CSC a cos /3 + ^ — = -^r, 



{l + d')jp[ 



(4.86) 



w — 



i\2i-iru{i+£)(}: 



sin a sin /3 



a'^' - {a" + (-l)V(l + d'')a'(5[)i^ 
^ 2da'{d cos o; sin /3 + cos /3) {y + dZ) cot a ^ Q;''(d cos a sin /3 — cos /3) ((i3^ — Z) 



1 + ^2 



(1 + c/2)q;' 



[{id? - l)y - 2dZ) COS a sin/3 + ((1 - c/^)^ - 2dy) cos/3]/3i' a' 
H ^ ;7rTT:i A sma sm p 

+q:'A'(cosq: sin /3 + (icos/3) — a!'3^(cos a sin/3 + (icosa cos/3) 



. ^ (a')2e(-i)''''(i+'^')^i(cosa sin/3 + (icos/3), 
-da'Z CSC a sm /3 + ^ —^r 



(1 + ^2)^/5; 
29 



(4.87) 



and p is given in (4-81). 

Case 2. 72 = ai = 0, (7T)' + - (-1)^7^ = and 71 ^ 0. 

According to (4.51) and (4.54), 

% = A;2 = 0, hy ^ -2(3' cosa. (4.88) 

7i 

Note 7 = 7^. Moreover, (4.52) says 

kx = 0, /iA^ = -2/3' sin a. (4.89) 

Furthermore, (4.61) yields 

fx^-. (4.90) 
71 

Besides, (4.40) implies 

hz^-{fx + ky) = Q. (4.91) 
Under the condition (4.60) and (4.61), (4.58) and (4.59) are equivalent to 

/^ = % = 7i = 0, -{-ly^v^ + i^' = (4.92) 

Write 7i = ai as a real constant. We have: 

a = ase^-^)''"'^?*, ^ = aie^-^)"""?*, r = a^'e^-^)"""?* (4.93) 

for 02 G M (cf. (4.47)). By (4.60) 



-a" 



that is. 



-a" 



Thus 

U = -a' - (3'Z sin a + aie(-^)'"'^"?*(i?^ + asJ^^^ - - aia2A^^) , (4-96) 

V = a'A'-/3'Zcosa + e(-^)'''"?*(C,- + 2aia20^A'), W = -/3'(A'sina + :Vcosa). (4.97) 
Moreover, 



a" 



Ri = -— + a^e(-^)''2''"i*[(4Wo,r- + 051,^)0^" + 3a2(25o,^(a6 + c) + (5i,^ccos 6)^'^ 
+2o^a2(4a5o,r + (^i,r)'^^ + fi2(2(5o,r(Q^ ~ c) + 5i^rCsin6)3^ + (^i^^cci sin(6 — h\) 
+2(5o,r(&Ci — 61 c) + aia2(25o,r(ci — a6i) — (Ji^^Ci sin 61)^"] + (a'/?' cos a — (/?' sin q;)')Z 
-{[plf + (/3'sinQ;)2)A' + e(-i)'''^'^i*(aia20 - 2a'))(Cr + 2aia2e(-i)'''^'^?Vr'^), (4.98) 
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a" 

-^X + (3(/30' cos^ a - [olf)y - ((/?' cos a)' + sin a)Z 
^axe^-^^''^'"''^\irC,r + 2aia2^r^r'^ + 2a\al<^ri}rX'^), (4.99) 

i?3 = {a cos q; — sin (x)')X — cos a)' + q;'/3' sin q;)3^ + a Z esc 2q;. (4.100) 
By (3.12), 

V = - a^e(-^)^^-^^[ ^^^'^°--+f^'-^'^' + a2(25o,.(a6 + c) + (5i,,ccos 6)^'^ 

2 2 

+a^a2(4a(^o,r + <^i,r)-;^ + a2(2(5o,r(af' ~ c) + ^i^^csin b)Xy + (5i,,.cci sin(b — 6i) 

aia2(2(5o,^(ci - a6i) - 5i,^ci sin 6i)A'2 ((a')^ + (/^'sina)^)^'^ 



+25o,.(6ci - 610))^" + -^--v--u,.v-. ] ^ ^ 

+ ((/?' sin a)' - a'ii' cos a)XZ - e^-^'^^ '"'^^\aia2<p - 2a'){irX + aiase^-^)'''^"?^,;!^^) 
_(_l)-e(-i)''^'^?*[a2e(-^)''-i^[aia2:^^0,V'r - 2-^02(01(40^0,1 + ^i,.)^'' + (fl) 

-\-'&rC,r - ((^0,l2(a&l + Ci) + 5l,rCi)3^] + a' + 02>'V'r " ar^<^2'?!>r)] ^ ^ 

(3(/3')^cos2a- (a')^)3^^ ^, a"Z2csc2Q;, 
V > +((/3^cosQ)^ + Qi^/3^sina)yZ } (4.101) 



2 

where modulo the transformation in (1.10), where 

£^ = h^e^ + cie~^ , £1 = ci cos(ro + 61) (4.102) 

in connection with (4.56) and w — oi3^. 

By (3.3) and (3.5), we have the following theorem: 

Theorem 4.3. Let a he a function of t and let 0,01,02,6,61,0,01 he real constants. 
Denote {3 as in (4-91). Define the moving frame X, y and Z hy (3.1) and (3.5), and 

00 = e"i^ - oe-"i^, 01 = sin(oi:V), -00 = e^^^ + oe""!^, -01 = cos(oiJ^), (4.103) 

^0 = ^e-^i^ - ce""^^, 6 = csin(oi3^ + 6), Co = ^e"^'^ + ce"'^^^, (4.104) 
Ci = ccos(oiJ^ + 6), i?o = 6ie"^^ - cie""!^, i?i = ci sin(oiJ^ + 61), (4.105) 
£0 = 6ie"i^ + cie""!^, £1 = cicos(oi3^ + 6i). (4.106) 
For r = 0, 1; we have the following solution of Navier- Stokes equations (1.1)-(1.4): 

u = [-a'y + oie(-^)'''"i*(i?^ + 02:^0^ - CrX - oi02V'r'^^)] cos a 

-[a'X + e(-^)'"''"?*(C^ + 201020^^")] sin a, (4.107) 
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V = [-a'y + aie(-^)'"'^"?*(i?^ + a23^(/>^ - Cr^ - aia2AX^)] sin a cos /? - P'Z cos P 
-[a'X + e^~^)'"''"i*(^^ + 2aia2(j)rX)] cos a cos /? + sin a + cos a) sin (4.108) 

w = [-a':^ + aie(-^)'''^"?*(i?^ + OaJ^^,. - - aia2A'^^)] sin a sin /3 - /3'Z sin (3 

- [a'X + e^-^^^'^'^^^Cr + 2aia2(l)rX)] cos a sin /? - sma + y cos a) cos /? (4.109) 

and p is given in (4-101). 
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